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Abstract 

The crystalline ferromagnet with dipole-dipole interactions is studied within the framework of the three- 
dimensional random-field Ising model. The field distribution function and the analytic expression for temperature 
dependent free energy are found. The dipolar induced magnetic anisotropy energy of bulk cubic ferromagnetic 
materials is calculated. 

In recent years several papers describing the magnetic anisotropy induced by the dipole-dipole interactions in 
ultrathin films were published In Refs. ||,|| the system consisting of several ferromagnetic monolayers was 

considered in the framework of a classic Heisenberg model. The temperature driven transition from out-of-plane 
to in-plane magnetization was explained. In Ref. Q the magnetic anisotropy of the two-dimensional lattice was 
analyzed by means of the dipole-dipole interactions. 

The random-field method in the Ising model [Q allows us to add the dipolar interaction and to estimate induced 
magnetic anisotropy. 



1 Model 



Consider N + 1 spins positioned at the sites of crystalline ferromagnet lattice. Let us suppose that one of the site 
is located at the origin of the coordinates. The probability of finding the interaction field H to lie in the interval 
(H, H + <2H) at the origin is (§§ 
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where h^, (mj., r^) is interaction field induced by the atom with the magnetic moment m^. situated at the site r^, 5 
is Dirac's delta function. Then the random field distribution function will be 
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where r k (m k ) is the distribution over atom magnetic moments m k function. In the approach of the Ising model, 
when the magnetic moments are aligned with the selected direction given by unite vector s and their magnitudes 
are equal to too, r k takes the following form 

Tk (m fe ) = a k S (m s - m fc ) + (3 k S (m s + m fe ) , 

where a k and f3 k are the probabilities of finding magnetic moment to be oriented along and opposite to chosen 
direction. Let us consider a case where the interaction field h k (rrifc,rfe) is a sum of exchange and dipole-dipole 
interaction: 



u / \ T f \ mfc , 3(r fe • m fc )r fe - r^m fc 
h k (m fc ,r fc ) = J(r k ) + -i '— * , 

m k r° k 

where J(r k ) — J being exchange constant if kth atom is nearest neighbor and J(r k ) = otherwise. 
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Let us replace a k and (3k by their average values a and (3 which are given by the equation 

exp (mgs • H/fesT) 



2cosh(rn s • H/fc B T) 



W(U)dH, a + (3=1, 



(3) 



where fc^ is Boltzmann constant, T is temperature. 

The characteristic function A(p) = J W(H) exp(ip ■ H)dH, being the Fourier transform of the function 1T(H), 
takes the form 



A(p) = Y[ [aexp(ip ■ h fe (m s,r fe )) + (3 exp (-ip ■ h k (m Q s, r k ))] 



(4) 



After the inverse transformation the distribution function is 

w(H) = £ c i E c i 2 ■ ■ ■ E ctdF^p"-^^ 
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where C{ k = 1 are binomial coefficients. 

One can obtain the free energy per atom using relation F = — k B T hi J exp ( "jj. ^* ) W(H)dH. After integrating 



F will be 



F = -fcsT^ln 
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and substituting (||) into Eq. (||) F takes the form 



F = —k B T^^ In < a exp 
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Note that the terms describing contribution of the external field H cxt and single-ion anisotropy — mos • H cxt + 
D(sx + s* + s D being single-ion cubic anisotropy constant, may be added to the Eq. (Q). One can calculate a 
and /3 using equation derived from (||) taking account of the exchange coupling only: 



= 2a - 1 = Cz® j ( l ~ a) Z ~ 3 tanh 
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where m = a — /?, Z is number of nearest neighbors. 



2 Results 

Let us treat dipole interaction as small parameter with respect to exchange coupling (toq <C J), also consider the 
case of high temperatures (mg <C fcgT). Expand F into a series in powers of m\ and direction cosines s x , s y , s z . 
Taking the identity for cubic crystals 
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into account and retaining first nonzero anisotropic terms of this expansion we will obtain the anisotropic part of 
free energy 



F m = -36a(3 
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a being lattice constant, 



^= E 

j ,k,l — — OG 
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can be calculated numerically. For bcc lattice it was calculated that L 1 « 0.00454, i 2 « 0.04343, L 3 = ^ and for 
fee lattice Li » 0.03200, L 2 « 0.12912, L 3 = |. 

As one can see at T = ferromagnet is totally ordered and F an derived from Eq. (Q) is equal to zero since 
quadrupole moment of this system is zero. With increase in temperature l-Fknl raises and at T > T c F an ^ 
because of used approximation of Ising model (magnetic moments of all atoms are parallel to chosen direction) . 

The method proposed can be used for estimation of the magnetic anisotropy induced by dipole-dipole interactions 
in two-dimensional and bulk materials. 
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